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A Chilean seismic regionalization through a Kohonen neural network 
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A study of seismic regionalization for central Chile based on a neural network is presented. A 
scenario with six seismic regions is obtained, independently of the size of the neighborhood or the 
reach of the correlation between the cells of the grid. The high correlation between the spatial 
distribution of the seismic zones and geographical data confirm our election of the training vectors 
of the neural network. 
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I. INTRODUCTION 

The idea of seismic regionalization can be traced back 
to 1941 when Gorshkov [l[ published one of the first 
studies on seismic regionalization for the URSS. Based 
on these earlier studies, Richter 0], in a paper entitled 
"Seismic regionalization" , founded the basis of a method- 
ology that allows to have information dealing with pre- 
venting the potential damage that an earthquake may 
cause through the definition of hazard zones. 

In a seismic country like Chile, and based on the ur- 
gent goals established by Richter, a proper regionaliza- 
tion map is desirable. Attempts to fill the gap in these 
matters have been done from time to time. In [3| is found 
probably the first regionalization study for Chile where 
the authors Gajardo and Lomnitz used the same methods 
of seismic correlation Tsuboi used in the case of Japan. 
Later works on the subject were mainly focused on en- 
gineering aspects as Walkner Labbe [f| and Barri- 
entos [6|, discussing on seismic risk. Essentially, these 
studies used the technique proposed by Cornell [3] and 
Algermissen and Perkings [8|. A more recent study was 
made by Martin [9[. Here the author used a longitudi- 
nal division criteria for the country using the line that 
separates the deep and superficial earthquakes, which 
leads to a high degree of coupling between the Nazca 
and Southamerican plate. This criterium enable him to 
identify two macro-zones in the continental zone: In the 
coast, with hypocenters with depths lesser than 40km, 
and the mountain macro-zone. As a criteria for latitu- 
dinal division, Martin computed the b parameter of the 
Gutenberg- Richter law for an initial surface, and then he 
varied the surface iteratively. If b changes dramatically, 
it means that a seismic zone has been crossed. These are 
the nine regions Martin found: in the coast; three zones. 
In the deep mountain zone, he distinguished four zones. 
In the surface mountain zone, he distinguished a single 
zone. Finally, the method determined the aseismic zone 
of Magallanes. 

In this work we implement a Kohonen neural network 
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(from here on NN) to determine the seismic zones. In 
general, a neural network [l2j can be defined as a com- 
puting system made up of a number of simple, highly 
interconnected processing elements, which processes in- 
formation by their dynamic state response to external 
inputs. NNs are typically organized in layers. Layers are 
made up of a number of interconnected 'nodes' which 
contain an 'activation function'. Patterns are presented 
to the network via the 'input layer', which communicates 
one or more 'hidden layers' where the actual processing 
is done through a system of weighted 'connections'. The 
hidden layers then link to an 'output layer' where the 
answer is output. Most NNs contain certain types of 
'learning rule' which modify the weights of the connec- 
tions according to the input patterns which are exposed 
to. 



In this work we use a special NN appropriated to dis- 
criminate spatial distribution. A Kohonen NN [l8j](or 
as it is usually called a self organising map) is an arti- 
ficial neural network that is trained using unsupervised 
learning to produce a low-dimensional (typically two di- 
mensional), discretized representation of the input space 
of the training samples, called a map. The network must 
be fed a large number of example vectors that represent, 
as close as possible, the kinds of vectors expected dur- 
ing mapping. The examples are usually provided several 
times. In this way, the NN finds by itself the best distri- 
bution of zones, based on the training vectors, without 
any other specification than geographical characteristics. 

In the context of seismic hazard or seismic risk, there 
have been many attempts to use this technique. For ex- 
ample, in [l5| the authors used artificial NN to discrim- 
inate between earthquakes and und erg round nuclear ex- 
plosions. Also for seismic detection [16[, in which the NN 
is trained to reco gniz e signal patterns, and also for earth- 
quake prediction [13]. This work is organized as follows. 
The next section discusses the election of the training vec- 
tors or training set, from which the NN starts its learn- 
ing process. Then, in section III we discusses the way we 
perform the training process. Section IV is devoted to 
the analysis of the results and we finish the paper with 
conclusions. 
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II. TRAINING SET 

The election of the training set is the crucial step to- 
wards the implementation of a useful NN. The usefulness 
of a NN depends on what we are interested to get from 
it. As we have discussed in the previous section, the ne- 
cessity is mainly the one of a hazard categorizer in the 
Chilean territory. In this section we discuss the criteria 
in the election of our training set. 



A. Database 

The idea was to associate the earthquake information 
by zone, frequency and magnitude. We expected the NN 
finds by itself what we call here 'zone', by introducing in 
it specific information about its seismic history. In this 
work we use as a database all the earthquakes classified 
with Richter magnitude equal or larger than 4.5 Ms in 
the interval 1957 to 2007. The sources for the data we 
used, were obtained from the catalog CERESIS [l3j and 
from the database of the USGS (from 1973 to the 
present) . 

The area covered by our data set is Chile's continental 
territory included inside a rectangle defined around the 
meridian 72° W as the North-South axis, between the 
parallels 17° S and 56° S (see the figures for details). 

The interval of fifty years ensures us to have a very 
representative set of data. It is long enough to have 
earthquake frequency histograms with representative in- 
formation which could be projected to a larger period of 
time. 



B. The Kohonen training vectors 

In order to classify the zones, we divided the area se- 
lected in the previous section by 156 rectangular cells of 
1° x 1°. Each cell was characterized by a seven dimen- 
sional vector, indicating the observed historical seismicity 
between 1957 to 2007. The 7D vector is explained below 
with more details: 

x\\ Mean deep of the seismic sources 

X2- Number of earthquakes with magnitude > 4.5 Ms 

X3: Number of earthquakes with magnitude > 5.5 Ms 

X4: Number of earthquakes with magnitude > 6.5 Ms 

£5: Maximum observed magnitude 

x§: Central horizontal coordinate of the cell 

£7: Central vertical coordinate of the cell 

These are considered the minimum number of sensitive 
parameters which characterize each event. No geophysi- 
cal input was necessary. 
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FIG. 1: Meta stable configuration after 1000 iterations of the 
neural network. This result was obtained independent of how 
many classes we ask for, revealing a strong . 



III. TRAINING PROCEDURE 

The 156 vectors were introduced into the Kohonen NN. 
The free parameters in the running are two: the number 
of iterations and the number of classes or regions. By 
inspection, it was observed that in all the runnings, and 
after an average of 2500 iterations, the NN reached a sta- 
ble configuration, indicating that the system had found 
a local minimum. 

In an initial training session, we forced the NN to cre- 
ate 5, 10 and 15 classes. In all these cases, it was observed 
that during the first 1000 iterations, the system reached a 
stable configuration with just three classes, as it is shown 
in Fig. |1]). 

We called these three classes "the fundamental 
classes" . Let us explain in detail each one of these. Ac- 
cording to Fig. |T]), in the northern part of the coun- 
try, we distinguish two classes or two seismic zones of 
high seismicity, associated to the interaction between the 
Nazca and South-American plate. In the south, we dis- 
tinguish another class or seismic zone, with low seismic- 
ity, associated to the interaction among the Southameri- 
can, Scotia and Antartic plates. 

Enhancing the detailed study, we forced the NN to cat- 
alog six classes or seismic regions. In this case, after 2500 
iterations, the final result was the same, essentially inde- 
pendent of the initial size of the neighborhood. Figure 
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FIG. 2: Transient stable configuration after 2500 iterations 
and independent of the initial reach of the correlation (or 
initial size of the neighborhood). The numbers are explained 
in Table I. 



TABLE I: This is a summary of the results for the parameters 
in the Gutenberg-Richter law. The equation reads log N = 
a + bM, where M is the magnitude in Ms units. It also is 
given the x 2 01 the fit and the maximum magnitude (in Ms 
units) for each region. 



Region 


a 


b 


x 2 


Max.Mag. (Ms) 


N(> 4.5) 


1 


3.844 


-1.069 


0.9073 


5.5 


0.11 


2 


4.209 


-0.922 


0.8650 


6.9 


1.15 


3 


3.298 


-0.699 


0.9373 


7.7 


1.42 


4 


2.702 


-0.520 


0.9721 


8.3 


2.30 


5 


3.784 


-0.685 


0.9405 


7.1 


5.02 


6 


3.702 


-0.555 


0.9388 


7.8 


16.03 



TABLE II: Here it is displayed the seismic zones under the 
hierarchy of a growing order according to the |6| parameter 
of the Gutenberg-Richter Law. NN used as input data the 
amount of seisms observed in the last fifty years, we assume 
that the final hierarchy must be mainly based on the num- 
ber of earthquakes which are equal or larger than 4.5 Ms, 
according to Table I. 



Hierarchy 



Zone 



\b\ 



0.520 
0.555 
0.685 
0.699 
0.922 
1.069 



associated to each region depends on the chosen vectors 
representing each cell. 



B. Gutenberg-Richter Law 



([2|) shows the final result using 6 clusters or categories, 
2500 iterations and an initial size of the neighborhood of 
5. 



IV. ANALYSIS 



For each region we found the parameters of the 
Gutenberg-Richter (GR) law. Here A^(> M) is the mean 
amount of earthquakes of magnitude larger than M ob- 
served in a year, in a normalized area (equivalent to a 
cell of 1° x 1° in the equator). 



A. Island regions 

As we can see in Fig.©, in certain regions some dis- 
continuities occurred, through an interruption of another 
region. We called these island regions. 

In the north part, region 6 becomes separated by two 
subregions; region 6A (coast) and region 6B (mountain). 

The north part, region 5 becomes separated in three 
subregions; region 5 A (including Iquique city), region 5B 
(including Copiapo city) and region 5C (centered in San- 
tiago city). 

Although they are discontinuous regions, the Kohonen 
NN associates these in the same category because they 
presented similar seismicity. For that reason we stress to 
name it just region 5 and 6. Of course, the seismicity 



V. CONCLUSIONS 

Considering the semi-logarithmic GR law slope, we ob- 
tain the hierarchy of the six regions, as it is shown in 
Table II. 

In Table I, and looking at the first and last columns, we 
observe the hierarchy based on the seismicity of the zones. 
As it was explained in section III, the NN was trained 
with information about the number of earthquakes per 
range. Then, the definition of seismicity deduced from 
the NN must be associated to the magnitude and angular 
orientation of the 7D vectors, which strongly depends on 
the number of earthquakes equal or larger than 4.5 Ms. 

In this context, zone 1 is an aseismic zone, and zone 6 
is the most seismic one. 
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It is possible that zones 5A and 6B are parts of a larger 
one extending themselves towards Peru and Bolivia. 

The zones 5B and 5C have an interruption between 
latitudes 29° and 32°, which coincides with a geographic 
accident; the intermediate depression despair, being re- 
placed by the fusion between the Andes and coast moun- 
tains. Then, zone 5 is strongly associated to the inter- 
mediate depression and zone 3 is strongly related to the 
Andes mountains. 




FIG. 3: North and South part of Chile. Here is the division 
of discussed seismic zones in the text. 



In the north part of the country, the high seismicity 
observed can be understood as the result of the activity 
of the subduction between the Nazca and Southamerican 
plates (at a rate of convergence of ~ 10cm per year, as it 
is well known). 

In the south, we obtain a low seismicity in agree- 
ment with the well known result of a convergence of the 
Southamerican, Scotia and Antartic plates, at a rate of 
~ 2 cm per year. 

In this work we have demonstrated that using a Koho- 
nen NN, and the database available up to date, it can be 
possible to differentiate six seismic zones for continental 
Chilean territory Some of them are apparently discon- 
nected zones (we called these here island zones) , but this 
is in fact an effect produced by considering just an ar- 
bitrary portion of territory (that of the political Chilean 
one) in the analysis. Also, the division between these 
island zones can be clearly associated with major geo- 
graphical accidents (as the intermediate depression and 
some deep valleys) that were not used as input data in the 
NN. All these facts show the power of using NN analysis 
in seismic regionalization. 

Acknowledgments 

JR wants to thank TGT for the support through grant 
number 2122. VHC wants to thank Rafael Valdivia for 
useful discussions. 



[1] G.P. Gorshkov, Publi. Bureau Central Seismologique In- 
ternational. A19, 25 (1956). 

[2] C.F. Richter, BSSA, 49, 123 (1959). 

[3] E. Gajardo and C. Lomnitz, Seismic province of Chile, 
Publication 11, 1529 (1958). 

[4] P. Welkner, Estudio de la Sismicidad en Chile y su apli- 
cacion al calculo antisismico, Thesis, Geophysics Depart- 
ment, University of Chile (1964). 

[5] J.C. Labbe, Relaciones macro sismicas para la evaluacion 
del riesgo sismico en Chile y California, (1976). 

[6] S.E. Barrientos, Regionalizacion Sismica de Chile, M. Sc. 
Thesis, School of Engineering, University of Chile, 1980. 

[7] C. Cornell, Seism. Soc. Am. Bull, 58, 15831606 (1968). 

[8] ST. Algermissen and D.M. Perkins, Geolog. Surv. Op.- 
File Rep. 76-416, (1976). 

[9] A. Martin, Hacia una regionalizacion y calculo del peligro 
sismico de Chile, Engineering Thesis, School of Engineer- 
ing, University of Chile, 1990. 
[10] G. Giacinto, R. Paolucci and F. Roli, Pattern Recognition 



Letters, 18, 1353 (1997). 
[11] M. Sanchez-Silva and L. Garcia, Earthquake Spectra, 17, 
89 (2001). 

[12] R. Rojas, Neural Networks - A Systematic Introduction, 

Springer- Verlag, Berlin, New- York, 1996. 
[13] CERESIS Catalog. Catalogo de terremotos para America 

del Sur (1985). 
[14] USGS. Web site: http//earthquake. usgs.gov/ 
[15] F.U. Dowla, SR. Taylor and R.W. Anderson, BSSA, 80, 

1346 (1990). 

[16] J. Wang and T.-L. Teng, BSSA, 85, 308 (1995); W.- 
Y. Kim, V. Aharonian, A. L. Lerner-Lam, and P. G. 
Richards, BSSA, 87, 569 (1997). 

[17] G. Giacinto, R. Paolucci and F. Roli, Pattern Recogni- 
tion Letters, 18, 1353 (2000). 

[18] T. Kohonen, IEEE Transactions on Computers, Vol. C- 
21, 353 (1972); T. Kohonen, Self-Organization and As- 
sociative Memory, Springer- Verlag, Berlin (1984). 



